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Abstract 

We study the impact of the two-loop electroweak corrections on the production of a 
Higgs boson via gluon-fusion in proton-proton collisions at LHC energies. We discuss 
the prescritpion to include the corrections to the hard scattering matrix element in 
the calculation of the hadronic cross-section a(p + p — ► H + X). Under the hypothesis 
of factorization of the electroweak corrections with respect to the dominant soft and 
collinear QCD radiation, we observe an increase of the total cross-section from 4 to 8%, 
for m H < 160 GeV. This increase is comparable with the present QCD uncertainties 
originating from hard scattering matrix elements. 
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1 Introduction 



The Higgs boson is one of the missing ingredients of the Standard Model and its discovery 
represents one of the most important physics goals of the LHC. This goal will be achieved 
only if we can predict with high accuracy all the production cross sections of this particle 
and if we understand in detail the different decay channels and the relative backgrounds. 

At the LHC, the gluon- fusion is the dominant production mode over the entire range of 
interesting values of the mass of the Higgs particle (lOOGeV < m# < ITeV). In particular, 
in the range lOOGeV < m# < 2m t this production mode is larger by almost one order 
of magnitude with respect to the next important mechanism, the vector boson fusion. It 
is, therefore, very important to have a precise prediction of its cross section and a reliable 
estimate of the remaining theoretical accuracy. 

The total cross section for the Higgs boson production by gluon fusion in the LO ap- 
proximation was calculated in the late seventies [T]. It is an 0(a$a) calculation, since the 
Higgs couples to the gluons only via a heavy-quark loop (the most important contribution 
is the one due to the loop of top). For what concerns the higher orders, the calculation of 
the NLO QCD corrections have been done in the infinite m t (mass of the top) limit in [2], 
and, with the full quark mass dependence, in [3j. Besides of the fact that the infinite m t 
approximation should be valid in the Higgs mass range m H < 300 GeV, it has been noticed 
[3| that this approximation works also for values of mg beyond the top threshold, and up 
to masses of O(lTeV). The total effect of the NLO QCD corrections is the increase of the 
LO cross section by a factor 1.5-1.7, giving a residual renormalization/factorization scale 
dependence of about 30%. The electroweak corrections were evaluated in the infinite m t 
limit in [H] and turned out to amount to less that 1%. The attention was driven by the 
evaluation of the NNLO QCD corrections, carried out in the infinite mt limit by several 
groups pa|. The calculation shows a good convergence of the perturbative series: while the 
NNLO corrections are sizable, they are, nevertheless, smaller that the NLO ones. More- 
over, the NNLO corrections improve the stability agaist renormalization/factorization scale 
variations. The effect due to the resummation of soft-gluon radiation has been included in 
[Zj, and the remaining theoretical uncertainty, due to higher-order QCD corrections, has 
been estimated to be smaller than 10%. Finally, several efforts were also devoted to the 
calculation of QCD radiative corrections to less inclusive quantities, such as the rapidity 
distribution, recently evaluated at the NNLO |Hj, or the transverse momentum (g^) distri- 
bution P|, which, in [THI, is evaluated using the fixed-order perturbative results up to NLO 
in QCD and the resummation up to the NNLL. 

Motivated by this accurate scenario, the NLO electroweak corrections to the gluon fusion 
were reexamined recently. In [TTJ the contribution to the partonic cross section due to the 
light fermions were calculated. It turned out that they are sizeable. In particular, in the 
intermediate Higgs mass range, from 114 GeV up the the 2mw threshold, these corrections 
increase the LO partonic cross section by an amount of 4-9%. For larger values of the 
mass of the Higgs, ran > 2mw, they change sign and reduce the LO cross section; however, 
in this region the light-fermion corrections are quite small, reaching at most a —2%. In 
[T2] . also the remaining electroweak corrections due to the top quark were calculated as a 
Taylor expansion in m#/(4m^/). They are valid for m H < 2mw, range in which they have 
opposite sign with respect to the light-fermion corrections. However, the corrections due 
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Figure 1: Lowest order (a) and generic NLO-EW (b), (c), (d) Feynman diagrams. The 
solid lines are fermions. The wavy lines are gauge bosons (V = W,Z). 

to the top quark are smaller in size, reaching at most a 15% of the light-quark ones. 

The impact of the NLO electroweak corrections on the hadronic cross section has not 
been discussed yet. We present here the effect of their inclusion in the calculation at the 
hadronic level. 

2 Inclusion of the Two-Loop Electroweak Corrections 

The partonic gluon fusion process occurs, in lowest order, via one-loop diagrams, as the 
one depicted in Fig. 0(a); in the loop run only the top and the bottom quarks, because of 
the Yukawa suppression of the lighter quarks. The NLO-EW corrections are schematically 
represented by the diagrams in Figs. 0(b), (c) and (d). In particular, in Figs. 0(b) and (c) 
the WWH/ZZH couplings avoid the Yukawa suppression, and, therefore, in these diagrams 
the fermionic line represents all the possible flavours: light flavours, evaluated in [TT], and 
top quark, evaluated in [12]. In Fig. (d), instead, the fermionic line can represent only 
the top quark [T2| . 

At the hadronic level, we consider the Higgs boson production at the LHC, and therefore 
in proton-proton collisions. The hadronic cross section can be written as: 




= *o (l + K% CDonl y(a s (» 2 ),v 2 ,M 2 )) (1) 

where the partonic processes initiated by partons (a, b) are convoluted with the correspond- 
ing parton densities fi tP (x, M 2 ), {i = a, 6), evaluated at a scale M. The effect of the higher 
order QCD and EW corrections is described by the two functions K QCD ~ only and Sew, 
obtained by factorizing the lowest order cross section cr . 
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In the partonic cross section, QCD and EW corrections have been factorized. This 
ansatz is valid up to subleading higher order corrections which start at the 3-loop level 
(i.e. 0(aa s ) with respect to the lowest order). The factorization of the QCD initial state 
collinear divergences holds for the hard process described by the electroweak NLO correc- 
tions, following from general arguments of the factorization theorems and from the universal 
nature of the initial state collinear radiation. In fact, the whole set of EW corrections is 
characterized by a scale m w , much harder than the one typical of the leading collinear 
emission. In addition, in the limit of light Higgs, the EW corrections can be expanded as 
a Taylor series in powers of m H /m w and the EW corrections vertex becomes effectively 
pointlike. In this regime the factorization of the QCD collinear divergences becomes rigor- 
ous. For heavier Higgs masses, the factorization should still be valid only in leading order, 
due to the modifications induced by the EW form factor. 

3 Numerical Results 

The hadronic proton-proton cross section has been calculated at LHC energy, i.e. \/~S = 14 
TeV, in NNLO-QCD accuracy, i.e. setting 5 EW = 0, using the MRST2002 NNLO parton 
distribution functions [1 3j . The theoretical uncertainty due to the choice of the renormal- 
ization scale \i and of the factorization scale M has been canonically estimated by setting 
M — ii equal to m H /2,m H ,2m H respectively. The predictions, shown in Fig. [2] (dotted 
lines), vary by approximately ±8% with respect to the central value. This uncertainty is 
further reduced when including the effect of the resummation of all the initial state soft 
gluon radiation 0. 

The two-loop electroweak corrections have been added according to Eq. (JTJ and setting 
M = /i = m H . The light fermion corrections can be evaluated for any choice of m H , whereas 
the top quark contribution has been computed by means of a Taylor expansion and is 
limited to the region m H < 160 GeV. The hadronic cross section increases from 4 to 8%, for 
m H < 160 GeV. As we can observe in Fig. the effect of the electroweak corrections is an 
increase of the cross section by an amount which is of the same order of magnitude of the 
NNLO-QCD theoretical uncertainty, and possibly larger than the uncertainty estimated 
after the resummation of soft gluon radiation. The main source of uncertainty on the 
hadronic cross section remains in the accurate determination of the parton distribution 
functions of the proton. 

The effect of the NLO-EW corrections is of great interest, because it enhances the most 
important Higgs production mechanism and, in turn, affects the absolute number of events 
of all the Higgs decay modes. 

Following Eq. (JH), the NLO-EW corrections can be implemented as a simple rescaling 
of the QCD hadronic cross section. This multiplicative factor is presented in Table [HI as a 
function of m H and can be fitted, in the range 114 GeV < m H < 155 GeV, by the following 
simple formula: 

hw{m H ) = 0.00961 + 6.9904 • 10~ 5 m H + 2.31508 • 10~ 6 m\ . (2) 

The computation of the NLO-EW corrections to the gluon fusion process has been de- 
scribed in detail in [IH H2J. The analytical expression of the probability amplitude has 



3 




120 140 160 180 200 220 240 260 280 300 

m H (GeV) 



Figure 2: The cross section a p+p ->H+x, in pb, is plotted as a function of the mass of 
the Higgs boson, between 114 GeV and 300 GeV. The dotted lines describe the band of 
NNLO-QCD uncertainty, for three values of the QCD factorization/renormalization scale 
[A = m H /2,mH,2'm H . The solid line is the NNLO-QCD (// = m H ) with the two-loop EW 
corrections, according to Eq. |H). The two-loop EW corrections include also the top-quark 
effect, for m H < 155GeV, but only the light quarks contribution for larger values of m H . 
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Table 1: Rescaling factor 5ew as a function of the Higgs boson mass. 



been expressed in terms of Generalized Harmonic PolyLogarithms (GHPL) |14j and has 
been implemented in a FORTRAN routine 1 . The GHPL can be evaluated numerically in 
several different ways: by direct numerical integration of the basic functions, by power 
expansions or by solving the associated differential equations. We have checked that these 
fully independent approaches agree. 



4 Conclusions 

In conclusion, the calculation of the QCD corrections to the production of a Higgs boson 
via gluon-fusion has reached a very high level of accuracy; the inclusion of the two-loop 
electroweak corrections, whose typical size for m H < 160 GeV is larger than 5% and then 
comparable or larger than the QCD uncertainty, is highly desirable. The main source of 
uncertainty on the hadronic cross section remains in the accurate determination of the 
parton distribution functions of the proton. 
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